home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Cream of the Crop 1
/
Cream of the Crop 1.iso
/
PROGRAM
/
NRPAS13.ARJ
/
CHSONE.DEM
< prev
next >
Wrap
Text File
|
1991-04-29
|
1KB
|
53 lines
PROGRAM d13r8(input,output);
(* driver for routine CHSONE *)
CONST
nbins=10;
npts=2000;
TYPE
barray = ARRAY [0..nbins] OF real;
VAR
glinext,glinextp : integer;
glma : ARRAY [1..55] OF real;
chsq,df,prob,x : real;
i,ibin,idum,j : integer;
bins,ebins : barray;
(*$I MODFILE.PAS *)
(*$I GAMMLN.PAS *)
(*$I GCF.PAS *)
(*$I GSER.PAS *)
(*$I GAMMQ.PAS *)
(*$I RAN3.PAS *)
(*$I EXPDEV.PAS *)
(*$I CHSONE.PAS *)
BEGIN
idum := -15;
FOR j := 1 to nbins DO BEGIN
bins[j] := 0.0
END;
FOR i := 1 to npts DO BEGIN
x := expdev(idum);
ibin := trunc(x*nbins/3.0)+1;
IF (ibin <= nbins) THEN
bins[ibin] := bins[ibin]+1.0
END;
FOR i := 1 to nbins DO BEGIN
ebins[i] := 3.0*npts/nbins*exp(-3.0*(i-0.5)/nbins)
END;
chsone(bins,ebins,nbins,-1,df,chsq,prob);
writeln('expected':15,'observed':15);
FOR i := 1 to nbins DO BEGIN
writeln(ebins[i]:14:2,bins[i]:15:2)
END;
writeln;
writeln('chi-squared:':19,chsq:10:4);
writeln('probability:':19,prob:10:4)
END.